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Abstract. Data from KEK on subthreshold p as well as on tt^ and 
production in proton-nucleus reactions are described at projectile energies 
between 3.5 and 12.0 GeV. We use a model which considers a hadron-nucleus 
reaction as an incoherent sum over collisions of the projectile with a varying 
number of target nucleons. It samples complete events and allows thus for 
the simultaneous consideration of all particle species measured. The overall 
reproduction of the data is quite satisfactory. It is shown that the contributions 
from the interaction of the projectile with groups of several target nucleons are 
decisive for the description of subthreshold production. Since the collective 
features of subthreshold production become especially significant far below the 
threshold, the results are extrapolated down to COSY energies. It is concluded 
that an p measurement at ANKE-COSY should be feasible, if the high background 
of other particles can be efficiently suppressed. 



FAGS numbers: 24.10.-i, 24.10.Lx, 25.40.-h 

1. Introduction 

Subthreshold particle production in a nuclear reaction is understood as production 
below the energy threshold of the considered process in a free nucleon-nucleon (NN) 
collision. It is thus a nuclear phenomenon which may be explained by rather different 
assumptions on the properties of nuclear matter and on the interaction dynamics. 
Hereby it is an open question to what extent subthreshold particle production is 
governed by properties of the nuclear ground state wavefunction and to what extent 
by the dynamical properties of nuclear matter, not reflected in the ground state 
description. The problem is far from a final solution at present and evidently requires 
a systematic study of high-momentum transfer processes, among them subthreshold 
particle production. 

At COSY a research programme is in progress which is devoted to the systematic 
investigation of subthreshold particle production. Subthreshold production of 
K+ mesons PP was one of the foundations for the design and building of the ANKE 
spectrometer which allows for the measurement of K+ meson production cross 
sections under the condition of an extremely high background of other particles. 
Recently, first results on K+ meson production have been published. Moreover, 

it was proposed to extend the research programme to the subthreshold production 
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of K~ mesons. This made a considerable extension of the ANKE spectrometer 
necessary, because additional equipment for the registration of negative particles had 
to be installed. 

Following this line, in this paper the question is examined if the investigation of 
subthreshold p production is manageable at COSY energies. The highest energy of 
2.8 GeV reached at COSY up to now is far below the threshold energy of 5.63 GeV for 
p production in NN interactions. This makes the measurement a real challenge due 
to the low cross section and due to the high background of other particles. From the 
physical point of view, however, deep subthreshold production far below the threshold 
becomes especially interesting, because the collective aspects of the phenomenon can 
be expected to become more and more dominating. Indeed, the created mass in case of 
p production (1.876 GeV) is essentially higher compared to K+ production (0.67 GeV), 
and the energy is 2.83 GeV below the NN threshold if a p production experiment would 
be carried out at 2.8 GeV. Such a difference is significantly larger than the value of 
0.58 GeV achieved so far at ANKE in the K+ production experiments at an energy of 
IGeV 0. 

First measurements of subthreshold p production in proton-nucleus collisions 
at higher energies had been carried out |3 |H1 El ^1 a few decades ago. Then 
investigations in nucleus-nucleus collisions ^1 El El E] and more recent studies 
of proton-nucleus reactions El E| followed. Most of the descriptions proposed 
so far are based on transport calculations [T^l^l^l^l^l^l^l^l?7l 
thermodynamical considerations j29ll3()imil32j or multi-particle interactions [HSj- All 
these approaches consider the p spectra without any relation to measurements of the 
other reaction channels. At KEK ^Hl, however, the spectra of tt^ and mesons 
were measured together with those of the antiprotons. The present approach is 
distinguished by considering simultaneously all these reaction channels. This allows for 
a more comprehensive determination of the model parameters and makes a prediction 
for COSY more profound. 

It should be stressed that the investigation of subthreshold production in 
connection with the other reaction channels is also of great importance for the 
understanding of this process. Fragmentation of the target residue is such a 
distinguished channel having a large cross section. In a participant-spectator picture 
there seems to be rather a weak connection between the interaction of the participants 
leading to particle production and the excitation of the residual nucleus during this 
process. At subthreshold energies, however, the competition between these two energy 
consuming processes may heavily influence the cross sections observed. In [35 this 
question was considered for proton-induced subthreshold production of K"*" mesons. 
There, the interplay of subthreshold production and the fragmentation of the residual 
nucleus was investigated in the Rossendorf collision (ROC) model, which allows the 
treatment of hadronic and nuclear reactions in a unified way. Here, we extend these 
considerations to the case of subthreshold p production. Both the fragmentation of 
the residual nucleus and the interaction of the nucleons participating in the scattering 
process are treated on the basis of analogous assumptions. Even more important, the 
phase-space of the complete final state including the reaction products of the projectile- 
participant interaction as well as the fragments of the decay of the spectator system 
is exactly calculated. To the best of our knowledge this feature seems to be unique to 
the ROC model. 

In the ROC model the nuclear residue becomes excited during the reaction 
due to the distortion of the nuclear structure by the separation of the participants 
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from the spectators and due to the passage of the reaction products through the 
spectator system. In this global way final-state interactions are taken into account 
without making special assumptions concerning re-absorption, re-scattering, self- 
energies, potentials etc. for the various particle types. The discussion of these 
properties is rather controversial. If one considers e.g. re-absorption of p then the 
treatment reaches from no explicit consideration jSHllSHl over a global factor 0.1 for 
p-|-Cu [20! to taking the mean free path into account [23 • For Si-|-Si the survival rates 
of antiprotons are according to |28l 1241 IT!?] in the region of several %, while in [201 
a reduction of about 50% for p-|-Au and Si-|-Au from an analysis of p production at 
higher energies [JTj has been reported. The authors JS] conclude from an analysis of 
the target-mass dependence of their data that re-absorption of p is much smaller than 
expected from the corresponding pN cross section. Thus, it is far from being clear to 
what degree antiprotons are absorbed in the nuclear medium. 

The plan of the paper is as follows. In section|21the main ingredients of the ROC 
model are explained, which is used for the calculations to be presented. Section O 
contains a comparison of theoretical and experimental results for particle production 
and fragmentation with special emphasis on subthreshold p production. A prediction 
of the antiproton cross section at COSY energies and a comparison with the cross 
sections for tt^ and K~ production is made in section 0] A summary is given in 
section |31 

2. The model 

The ROC model is implemented as a Monte Carlo generator which samples complete 
events for hadronic as well as nuclear reactions. It makes no detailed assumptions on 
the intra-nuclear development of the interaction process, but calculates instead the 
statistical weights of the possible final states. The dynamics of the reaction is taken 
into account in form of empirical functions which modify the population of the final 
states. This is in contrary to transport models, where the interaction of the projectile 
with target nucleons and the subsequent interactions of particles originating from 
primary collisions with further target nucleons are modelled. The ROC approach 
does not need any parameterisation of elementary cross sections. Its applicability 
is not restricted at higher energies by the growing number of unknown elementary 
cross sections. All dynamic information is gathered in a few parameters which are 
either constant or change smoothly with energy and/or target mass. The energy 
necessary for subthreshold production stems from interactions of the projectile with 
few-nucleon groups, also called clusters in the following, and from the Fermi motion 
of these clusters. 

The model was successfully tested for pp interactions up to ISR energies in p^l!^ . 
while nuclear reactions were considered in the papers |84l 1401 1411 02| ■ In the following 
the basic ideas of the ROC model relevant for this paper are summarised. 

The cross section of the interaction of an incoming hadron with a target nucleus 
{A, Z) consisting of A — Z neutrons and Z protons is calculated as an incoherent 
sum over contributions from a varying number a of nucleons (thereof z protons) 
participating in the interaction 



A 



min{a,Z) 






az 



(1) 



a—1 z—max{0,a—Z) 



Subthreshold antiproton production in proton-carbon reactions 



4 



Here, s = denotes the square of the centre-of-mass energy of the projectile-target 
system with P — {E, P) being the total four-momentum, aaz stands for the partial 
cross section of the interaction of the projectile with a cluster (a, z), and the quantities 
dVF(s; cxaz) describe the relative probabilities of the various final channels cxaz- 

The partial cross sections Gaz account for sequential collisions between projectile 
and a nucleons of the target. Such a sequence is treated in the ROC model as a 
collision between the incident particle and a cluster consisting of a nucleons capable 
of sharing their energy in close analogy with the virtual clusters of the cooperative 
model gnilllllSlEnillZlEHlElEni- The aaz are calculated using a modified version 
of the Monte Carlo code which is based on a probabilistic interpretation of the 
Glauber theory jSl]. We use the profile function 

A 

W_PA{ri)A^ri 

i=l 

of the considered nucleus, which depends on the nucleon density pAi'f'i) and the 
probability 

Pi = exp{-d'^Tr/aNN) 
for an interaction of the projectile and the ith target nucleon with di being the distance 
between the interacting particles. The nucleon density j^Hl 

PA{r) cx (1 + r^ilMf ~ e')/f + e^V/^]) exp(-rV/2) (2) 

of light nuclei A < 20 can be derived from a standard shell model wavefunction with 
rj = {A- 4)/6 , p = 6^(1 - 1/A) and / = 1.55 fm. Then the NN cross section gnn 
is adapted such that the integral of the profile function over the impact parameter b 
reproduces the total inelastic pA cross sections 

< = / d^TAib), 

which are approximately constant in the energy region under consideration. The same 
calculation yields also the partial cross sections aaz we are interested in (for further 
details see 

The relative probabilities of the various final channels cxaz in ^ 

dW{s; aaz) cx dL„(s; aaz)pA{PR)T^{aaz) (3) 

are given by the Lorentz- invariant phase-space factor ALn{s\aaz) multiplied by the 
square of the empirical reaction matrix element T^{a.az) responsible for the collision 
dynamics. The Fermi motion is implemented via the momentum distribution of the 
residual nucleus pa{Pr)^ which is made a function of the number of participants a. It 
is taken as a Gaussian having a width of 

aa = ^a{A-a)/b/{A-l) pp (4) 

in accordance with the independent particle model |54| with pp being the Fermi- limit 
of the nucleus considered. No special high-momentum component |35l I55L I57L 1581 
15^ is used in this paper. 

The Lorentz-invariant phase-space is defined as the integral over the four- 
momenta of the n primarily produced final particles with energy-momentum 
conservation taken into account 

n n 

dLn{s;aaz) ^l[d''p^S{p^, - m^^) 6\P - Y,P^)- (5) 

i=l i=l 
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Here, the four-momentum of the i-th particle is denoted by = {si,Pi) withp^ — raf. 

For numerical calculations the 5-function in equation ^ has to be removed 
by introducing a new set of 3n — 4 variables to replace the 3n three-momentum 
components. It is reasonable to choose a set of variables, which reflects the underlying 
physical picture of the interaction process. Using recursion jBU] equation ^ can be 
rewritten in the form (see Appendix) 

dW{s- aaz) (X —^pA{Pn)dWR{Ml)mp{Ml). (6) 

Here, the four-momenta of the nuclear residue Pn — {En, Pr) and of the participants 

Pp = P-Pr = {Ep,Pp), (7) 

yield the corresponding invariant masses according to = P^ and Mp = Pp. The 
integral over the Fermi motion Pp p a{Pr) / 'i.E p separates the phase-space of the np 
nuclear fragments 

AWp{Ml) = dMlAL^, {Ml)TUa^,) (8) 

from the phase-space 

AWp{Ml) = dL„, {Ml)Tl{aP,,) (9) 

of the np ~ n — np final particles arising from the participant system. In (jSJ and 
the matrix element 

T^K,)^T|(afjT|(arj (10) 

is split into two factors describing residue fragmentation into channel a.^^ and 
participant interaction resulting in channel ct^^, accordingly. There is, however, a 
strong kinematic link between participants and spectators, since invariant mass Mp 
of the participant system, invariant mass Mp of the residue and the relative kinetic 
energy ^/s — Mp — Mp of these two particle groups are connected by energy- momentum 
conservation. For particle production to proceed the invariant mass of the participants 
must exceed the corresponding threshold value M^p which in turn depends on the 
number of participants. The heavier the effective target is the more energy is available 
for particle production. Another way to reach the threshold value goes via the Fermi 
motion, since the participant mass is a function of the momentum vector of the nuclear 
residue. The excitation energy of the target residue comes usually into play via the 
spectral function (see e.g. [SSI) derived from electron scattering data. Unique to the 
ROC model, however, is the treatment of the spectator system in close analogy with 
the participant subsystem. The ROC model calculates the complete final state of both 
the participants and the spectators. Thus, the huge amount of final channels of the 
spectator fragmentation infiuences directly the final state of the participant system 
and vice versa. 

The term dWp{Mp) in 10, © describes the interaction of the incoming proton 
with the group of participating target nucleons. Such a projectile-cluster reaction 
is treated in complete analogy to a hadronic reaction. In a first step intermediate 
particle groups called fireballs (FBs) are produced, which decay into so-called primary 
particles. The primary particles define the channels for which the weig hts 
are calculated. Among them are resonances, which decay subsequently into stable 
hadrons. In case of the presence of several nucleons among the decay products of a 
FB they may form nuclear fragments with equal probability for all possible channels. 
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The dynamical input of the reaction is implemented by the empirical transition 
matrix element 

T'pK.) = T? T^, Tl , (11) 

which describes the interaction process T-^ resulting in the production of a varying 
number N of FBs (N > 2), the production of hadrons T^^ via the creation of quark- 
antiquark (qq) pairs, the invariant-mass distribution of the FBs T^^, the transverse 
and longitudinal momentum distribution of the FBs, and, finally, some factors Tj. 
necessary for the calculation of the statistical weights. The interaction is assumed to 
proceed via colour exchange leading to the removal of valence quarks or of gluons from 
the interacting hadrons. Additional up, down and strange quark pairs are created in 
the ratio 

u:d:s = l:l:Xs (12) 

with As = 0.15. They form the varying number of FBs, which subsequently decay into 
the final hadrons. The transverse momenta of the FBs are restricted by an exponential 
cut-off (longitudinal phase-space) according to 

N 

t2 = J] exp(-7Pt,i) (13) 
1=1 

with the mean = 2/7. Two leading FBs, the remnants of the incoming hadron and 
of the cluster, get in the mean larger longitudinal momenta than the central FBs by 
weighting the events with 

= {XiX2)^. (14) 

Here, the light-cone variables Xi — (i?i +Pz,i)/(ef, +Pz,b) and = (-E2^^z,2)/(-Ep — 
Pz,p) are used with the four-momentum of the projectile given by pi, = (e?,, pf,) and that 
of the participants by Q. Each FB is characterised by two parameters, a temperature 
&FB and a volume Vfb- The temperature determines the relative kinetic energy of 
the particles the FB decays into via 

N 

TW^fb) = Y[iMi/eFB)KiiMi/eFB), (15) 
1=1 

while the volume defines the interaction region and influences mainly the particle 
multiplicity via the statistical factor T^^ oc Vpg~^. In l|15|) Ki stands for the modified 
Bessel function. Final hadrons are built-up by random recombination of the available 
quarks during the decay of the FBs. This procedure ensures automatically the 
conservation of all internal quantum numbers. Resonances decay later on until the 
final state consisting of stable particles is reached. For a more detailed discussion of 
the hadronic matrix element the reader is referred to |H9| . In this paper we use the 
same set of parameters as in (SHI for the description of the interaction of the projectile 
with a single target nucleon. Most of the other terms in QJ, however, contain clusters 
consisting of several target nucleons. The basic parameters of the FBs emerging from 
such an interaction are fixed by scaling the volume with the number of cluster nucleons 
a and the temperature parameter with a^^^^ according to 

Vfb = Vj^B a and Qfb = QTb" a"'^^- (16) 

It remains to consider the target residue, the structure of which is 
strongly disturbed by the projectile-participant reaction and subsequent final-state 
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interactions. This leads to the excitation and decay of the spectator system, which 
is characterised by the two parameters temperature Qji and volume Vr in the same 
way as the FBs emerging from the projectile-participant interaction. The part of the 
matrix element responsible for the residue fragmentation 

is identical with the corresponding factors Hll|l applied to the hadronic FBs. In order 
to restrict the excitation energy transferred to the residue we use the asymptotic 
approximation of (|15|l for large mass Mfj and small temperature 

TUQr) = VMr/Qr expi-Mn/eR). (17) 

An impact parameter dependence is assumed for the temperature parameter. This 
seems to be reasonable, because a peripheral collision with only few participating 
nucleons should excite the nuclear residue much less than a central collision with 
many participants. As a first guess we use 

Ofl = [1 - exp (-a/a^i/3)] (18) 

with d as parameter, here fixed to a = 0.5, which determines how fast the maximal 
temperature G^''^ is reached with increasing number of participants a. 

All factors still necessary for a correct calculation of the relative weights of the 
various channels are collected in complete analogy with the corresponding factors for 
the FBs in the term 

T,2,(afj = g(«fj -^^3 l[{2a, + 1) 2m.. (19) 

It contains the spin degeneracy factors (2(7^ + 1), the volume Vr in which the particles 
are produced with Vr = 4:Tt{A— a)R^/3 determined by the radius parameter Rr. The 
quantity gi^az) is the degeneracy factor for groups of identical particles in the final 
state of the residue decay and prevents multiple counting of identical states. 

Temperature and volume determine as in the case of hadronic reactions the 
number of final particles and their relative energy. The main difference consists 
in the value of the temperature parameter, Q™'jf « 300 MeV for hadronic and 
Qmax ^ 10 MeV for nuclear systems. New particles can be produced in hadronic 
reactions, while for nuclear systems, due to the much lower temperature, the nucleons 
of the initial state are recombined into various fragments without producing new 
hadrons. The volume parameter defines the distance between the fragments where 
the strong interaction ceases to work. This volume goes into the calculation of the 
weight of the final state. Due to the long range of the Coulomb repulsion Coulomb 
energy is still stored in the system at those distances. This energy is calculated in 
the Wigner-Seiz-approximation (see [HI]). The random distribution of the Coulomb 
energy among the charged fragments yields the final momenta of the fragments at 
large distances. 

In the calculation we assume that the fragments emerging from the residue decay 
are stable against particle decay. The possibility of hot fragments cooling down by 
subsequent particle emission is also implemented in the model (see PH), but it turned 
out that a readjustment of the volume and/or the temperature parameter yields quite 
similar results in both cases. So we decided in favour of the easier approach. Moreover, 
for a light nucleus as ^^C this question is of minor importance. 
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In concluding this section it should be mentioned, that it is the calculation of the 
relative probabilities dW{s;aaz), equation (O, which makes the difference between 
the cooperative IH [tSl EHl EH EHl EHl EOI and the ROC model, although the 
notion of clusters is similar in both approaches. In the cooperative model the relative 
probabilities of the various channels are calculated in the pure phase-space limit and 
the excitation and decay of the target residue is completely neglected. In contrast, the 
ROC model is far from being a pure statistical approach, because it contains specific 
information on the properties of nuclei and the interaction of the projectile with few- 
nucleon groups (clusters) in the nucleus. The model includes strong assumptions 
on the amplitude of the projectile-cluster interaction which lead to the restriction of 
excitation energy and momentum transfer to clusters compared to pure phase-space. 
A further nontrivial assumption concerns the various possibilities for the dissipation 
of the transferred energy inside the hit clusters either as relative kinetic energy of 
the decay products or as newly produced particles. In extreme cases the whole 
energy can be completely accumulated for the production of new particles. The other 
extreme is the emission of very fast secondaries with energies outside the kinematical 
limits of a projectile-(single)nucleon interaction without producing new particles. 
Also the supposition that clusters behave themselves like hadrons with respect to 
quark statistics is of great importance for the relative weights of the various final 
reaction channels. All these assumptions form a new approach to the consideration 
of local excitations and may, therefore, give essential and definite information on the 
properties of nuclear matter (excitation probability, distribution of deposited energy 
and transferred momenta, decay modes of few-nucleon fireballs etc). 

3. Comparison with data 

In order to fix the parameters describing the Fermi motion and the decay of the residue 
we consider the energy spectra of fragments from the reaction of 2.1 GeV protons 
with carbon measured at 90° (see figure For the volume a radius Rr = 1.7 fm 
is taken. This corresponds to about 1/3 of normal nuclear density and is within 
the region of break-up densities used for example by In figure Q the heaviest 

fragment, C, arises mainly from the quasi- free interaction of the incoming proton 
with a single target neutron and may gain energy only from the Fermi motion. Thus, 
the spectrum is highly sensitive to the width CTq , equation (gj , of the Gaussian for the 
Fermi motion. On the other hand, the lightest fragments are decay products of the 
residue and their kinetic energy in the laboratory is a superposition of Fermi motion 
and relative kinetic energy of the residue fragments. A satisfactory description of the 
fragment spectra in figure ^ is achieved using O^""^ = 12MeV, equation (fT5)) . and 
PF = 320MeV/c, equation Q). No attempt was made to improve the description 
of the spectra by implementing an additional high-momentum component, although 
especially the high-momentum parts of the spectra of the heavier fragments are 
underestimated. The spectator-participant picture is an idealisation, and not only 
the nuclear wavefunction but also secondary interactions between spectators and 
participants might be responsible for the observed deviations. All the following 
calculations are carried out with the above fixed values for 8^°^ and pp- It should 
be stressed in this connection that temperature @r and volume parameter Rn of 
the nuclear residue are highly correlated, a fact already observed in [3^] for the FB 
parameters in hadronic reactions. 

In figure 12 the momentum spectra of tt* and mesons as well as of p are 




Figure 1. DifTerential cross section as function of the kinetic energy of fragments 
produced in the reaction of protons with carbon 1621 (dots) compared to ROC- 
model calculations (histograms). 




Figure 2. Momentum spectra of tt^, and p 1181 (symbols) compared to 
ROC-model calculations (histograms). Experimental and calculated results for 
7r+ and K~ mesons are multiplied by the factors indicated in the legend. 
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Figure 3. Momentum spectra of tt^, and j3 I7UI (symbols) compared to 
ROC- model calculations (histograms). Experimental and calculated results for 
7r+ and K~ mesons arc multiplied by the factors indicated in the legend. 



compared to the KEK data ^^l- The overaU agreement is quite satisfactory in view 
of the large region of projectile energies and the variety of ejectile species calculated 
with one fixed parameter set. Particle yields are influenced by the suppression factor 
Xs — 0.15, equation p2|l . of strange quarks and by the algorithm for creating the final 
hadrons from the quarks produced in the first stage of the interaction process. Hadrons 
are built up in each FB independently according to the rules of quark statistics |M) 
by randomly selecting sequences of g's and g's. A qq gives a meson, while baryons 
or antibaryons are formed from qqq or qqq. From a given sequence of quarks the 
different hadrons are formed according to the tables of the particle data group j(i5j . 
There is no parameter which directly determines the ratio between meson and baryon 
production as e.g. in the PYTHIA-LUND model ,66,, _67,^, 69 . Only an indirect 
influence via the temperature and the volume parameter is possible, which change the 
relative weights of the FBs in dependence on their invariant mass and final particle 
multiplicity, respectively. 

It is therefore quite remarkable that the general trend of the data is well 
reproduced by using the parameters determined in |H9j from the consideration of 
hadronic interactions. The moderate increase with energy of the cross sections for 
the light mesons, the steep increase of that for the antiprotons as well as the shift 
of the maximum in the p spectra towards higher momenta is well reproduced by the 
calculations. 

At the highest energies the ROC calculations obviously tend to underestimate 
the data. Therefore, we compare in figure|31as a kind of cross-check a similar data set 
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Table 1. ROC results for the relative contributions to the total p cross sections 
at selected energies in % from terms containing various numbers a = 1 ... 6 of 
participants as well as for the total cross sections crp and the cross sections cr!?' 
with no freezing of kinetic degrees of freedom (see text). 



Energy(GeV) 


1 


2 


3 


4 


5 


6 


o-p(mb) 


o-f{mb) 


2.8 


0.0 


0.9 


51.3 


37.5 


9.8 


0.4 


7.810-1" 


6.3 10-" 


3.5 


0.05 


20.9 


56.6 


19.3 


2.9 


0.2 


1.510-6 


7.9 lO-'^ 


5.0 


10.9 


30.5 


43.1 


13.2 


2.0 


0.3 


2.9 10-* 


4.310-" 


12.0 


63.3 


21.7 


11.9 


2.6 


0.4 


0.05 


0.10 


0.13 



from [201 with ROC-model calculations. Here, the spectra are measured at a similar 
angle but in a much wider momentum region compared to the KEK data. In this case 
the tendency of the ROC results goes in the opposite direction and overestimates the 
cross-sections for the heavier ejectiles and p. 

The calculated cross sections are sensitive to the assumed temperature of the 
residual nucleus (see figureEJ- It should be stressed that the temperature derived from 
the fragment spectra yields also a reasonable reproduction of the p cross sections. This 
point is of special importance since it demonstrates that the link between the realms 
of fragmentation and of hadron production made by the ROC approach is obviously 
correct. 

The authors ^HI interpret their p data by using the "first-chance NN collision 
model" from 1851 where the internal nuclcon momenta were extracted from backward 
proton production j55| as a superposition of two Gaussian distributions. In this way 
the momentum spectra and the incident-energy dependence could be successfully 
reproduced by adapting one normalisation parameter. For the similar case of 
subthreshold K+ production it has been, however, argued [2] that the contribution 
from the high-momentum component should be negligible. 

In the ROC model the "first-chance NN collision" corresponds to the term in Q 
describing the interaction of the projectile with a single target nucleon. In table ^ 
the relative contributions of the terms with a given number of participants to the 
total p production cross section are summarised. It can be seen that the interaction 
with a single nucleon yields the main contribution only at the highest energy. The 
lower the energy the more contribute the terms with several target nucleons. At the 
lowest energy, 2.8 GeV, considered here the main contributions to p production arise 
from interactions with three and four target nucleons, while quasi-free collisions are 
negligible. At the highest energy, 12 GeV, there are still considerable contributions 
from proton-cluster reactions, although quasi-free collisions became the dominating 
process. Thus, the interpretation of the data by the ROC model is quite contradictory 
to the assumptions of a "first-chance NN collision model" . From the viewpoint of the 
ROC model the key quantity for understanding subthreshold particle production is 
the number of participating nucleons. A direct experimental determination of this 
number for subthreshold p production is highly desirable as discussed in 42 for the 
case of production. 

A certain similarity can be, however, found between the ROC approach and the 
multi- nucleon mechanism of j27j . where "the incoming proton and the interacting 
nucleons in the target act as sources of pions that merge to produce a nucleon- 
antiproton pair" . In this approach the number of participating nucleons obviously 
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plays a similar role as in the ROC model. 

Another interesting aspect of subthreshold particle production is the observation 
that the formation of light nuclei in the final channel alongside the produced particle(s) 
has to be taken into account for achieving a good reproduction of the measured 
cross sections. The formation of light nuclei leads to the freezing of relative kinetic 
energy of the nucleons, which is then available for particle production. We are 
speaking here about light fragments emerging from clusters, not about the decay 
of the residual nucleus. This effect has been discussed in 0H| for pion, in jHJj for 
photon, in |50| for K+ production in nucleus-nucleus reactions on the basis of the 
cooperative model, and in 05] j within the ROC model, for subthreshold production 
in proton-nucleus interactions. In order to demonstrate the importance of this effect 
for subthreshold p production we have in tabled in addition to the "normal" ROC 
results (Tp also the cross sections a?^ with no freezing listed. At the lowest energy 
considered the difference amounts to about one order of magnitude. This is relatively 
small compared to the three orders of magnitude reported by 021 for the reaction 
i2C(85MeV/A) + ^^C ^ Tr+ +X. The energy of 2.8 GeV is still far above the total 
production threshold of about 2.2 GeV for the coherent production at the whole target 
nucleus. Contrary to the energy considered in 0^1 also the number of competing 
channels is still rather high. This diminishes the influence of the neglected channels. 
The slight increase of cr?^ compared to ap at the higher energies is due to the smoother 
energy dependence of the relative probability dW{s; a.az) of the considered channel 
and the smaller denominator in equation as a result of neglecting a large number 
of final channels. 

In 1281 1241 I^l28j calculations have been carried out in the framework of transport 
theory, where questions concerning baryon self-energies, elementary p production 
amplitudes, p potentials and p re-absorption are of importance for the reproduction 
of the experimental data. In the ROC model there is no equivalent of these intra- 
nuclear properties implemented so far. Nuclear properties enter the ROC calculations 
only via the scaling 1)16(1 of the temperature and the volume of the FBs and via the 
excitation of the target residue caused by secondary participant-spectator interactions. 
Nevertheless, the spectra of all particle species are comparably well reproduced. This 
result should be further explored by considering heavier target nuclei, where the 
problem of secondary interactions becomes more important. On the other hand, the 
experimental investigation of p production on the lightest nuclei as ^He and ^i?e is 
of special interest, since the model predicts in the energy range between 3 and 5 GeV 
just the 3- and 4-nucleon groups as the main source of p production. Secondary 
interactions of the produced particles and the influence of the excitation of the target 
residue are practically absent in this case. This makes the interpretation of such data 
more clear than in the case of heavier target nuclei. 

4. Antiproton cross section at COSY energies 

Figure^ shows a ROC-model estimate of the antiproton spectrum for a bombarding 
energy of 2.8 GeV averaged over a polar angle 8 < 10°. In order to emphasise the 
experimental difficulties of a measurement at such a deep subthreshold energy, the 
spectra of the other negative stable particles are depicted, too. In the region of the 
maximum in the spectra around 500 MeV/c about ten orders of magnitude more 
7r~ mesons than antiprotons must be expected what makes a possible measurement 
in addition to the small p cross section a difficult task. 
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Figure 4. Momentum spectra of negatively charged hadrons calculated with the 
ROC model (histograms). 



In figure |S1 the energy dependence of the total production cross section for p in 
p + ^^C interactions is plotted. In the energy region considered the cross section 
changes by more than seven orders of magnitude. The importance of the temperature 
of the residual nucleus at low energies is demonstrated by a calculation with the 
temperature arbitrarily decreased by a factor of two. This increases the cross section 
by about half an order of magnitude at the lowest energy. With increasing energy the 
influence of this parameter diminishes and becomes negligible above 10 GeV projectile 
energy. The Fermi motion is of similar importance for subthreshold p production. A 
decrease of pp, equation Q), by about 30% makes the p cross section by nearly one 
order of magnitude smaller. Again the influence becomes negligible at high energies. 
These results underline once more how important it is to use a model which allows the 
verification of these aspects by independent measurements as has been done in this 
paper. It should be stressed that the selected parameter range does not reflect the 
uncertainty of the considered parameters which can be rather precise determined from 
the fragmentation data [HI] (see figure^. Instead, it is the aim to demonstrate the 
sensitivity of the results to the inclusion of the nuclear residue into the considerations. 
All results in the subthreshold region from approaches neglecting the excitation of the 
nuclear residue should be, therefore, taken with caution. 

But also the present prediction contains of course considerable uncertainties. 
From our experience with existing experimental data, in particular those shown in 
figures 121 and O we know that deviations between data and calculations are usually 
less than a factor of two. Since any drastic change of the model feasibility is not seen 
if the energy decreases from 3.5 GeV to 2.8 GeV we may expect the same level of 
uncertainty at this lowest energy. Certainly, for the present case of a prediction rather 
far below the region where experimental data are available it cannot be excluded 
that some processes or factors not included in the ROC approach may increase their 
significance with the energy reduction. A possible example could be the coherent 
production at the whole target nucleus jT^J. To be on the sure side we, therefore, 
assume at least the order of magnitude of the predicted cross section to be correct. 
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p + p + X 




projectile energy (GeV) 

Figure 5. Cross sections of p production as function of incident energy between 
2.8 and 12.0 GeV. The p cross sections calculated for three parameter sets (see 
text) are connected by straight lines. 



At the ANKE spectrometer an effective detection of particles produced at angles 
^ 10° in the momentum range 0.2 — 1.0 GeV/c is possible. Assuming a luminosity of 
10^^cm~^s~^ available with a carbon strip target and a cross section of 1 • 10~®mb 
expected at 2.8 GeV one gets a counting rate at the level of five hundreds events per 
week of beam-time. This estimate shows the feasibility of studying p production even 
at such a low subthreshold energy as considered here, if the experimental equipment 
allows for an efficient suppression of the negative mesons. 

5. Summary 

Subthreshold particle production is a collective phenomenon which is far from being 
completely understood. Data on subthreshold particle production can be reproduced 
within the ROC model by considering the interaction of the projectile with fcw-nucleon 
systems in complete analogy to the interaction with a single nucleon, also with respect 
to high-momentum transfer processes. It is the simultaneous consideration of the 
data for all particle types measured, which distinguishes the present approach from 
the previous attempts to interpret solely the subthreshold p spectra from KEK. The 
reproduction of the data requires no strong distortion of antiprotons in nuclear matter, 
which might be expected due to the fact that the free pA^ cross section is larger than 
the corresponding ones for the other particle species considered. This finding needs 
further confirmation by measuring p spectra at lower momenta than the region covered 
by the KEK data where the pN cross section increases further. 

In this paper the feasibility of measuring subthreshold p production at COSY- 
ANKE is demonstrated. It is argued that new insight into the mechanism of 
subthreshold production could be gained from an experiment nearly 3 GeV below the 
NN threshold, where the ROC-model calculations predict increasing contributions to 
the cross section from few-nucleon groups. The reliability of the predictions is verified 
by considering together with p production also the results for other particles (tt^, K^) 
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Appendix 

The n particles of the phase-space ^ of the final state of the reaction are divided 
into two groups, the decay products of the residual nucleus and the particles emerging 
from the projectile-cluster interaction. In a first step the decay products are separated 
by introducing their invariant mass and four-momentum Pr via the following 
identities: 

1 = J dMl 5{Pl ~ Ml) (A.l) 

1 = J d^PHSHPR~J2P^^ (^-2) 



with the definition 

Tin 



PR=J2p^. (A.3) 

i=l 

Inserting (|A.1|I . . . (IA.3|I into jSJ the phase-space factor becomes 
dL„{s) = dMld^PR 5{Pl ~ Ml) 

Y[d^p.6{p^-m^)6''iPn-J2p^) 



np 7ip 



Y[d'p,6{p^-m^)S*iPp-J2p,). (A.4) 

1=1 i=l 

Taking into account that for on-shell particles the identity 

d'p,Sipj-mj)^ f ^ (A.5) 
J 2ei 

is valid and considering the definition jSJ of the n-particle phase-space equation (IA.4p 
can be written as 

dLn{s) = dMl^ dLnAMl) dL„,(P|), (A.6) 

with Up — n ~ Up being the number of particles emerging from the participant 
interaction. The invariant mass of this participant group is given by Mp = ^ Pp 
with the four-momentum Pp = P — Pp. The separation into two groups makes 
sense, because the distribution of the momentum Pp can be derived from the internal 
momentum distribution p{Pp) of the target nucleus. 

Inserting the matrix element 1)10(1 and the momentum distribution p{Pl) into 
equation l(A.6(l and taking the definitions ^ and into account we arrive at the 
expression (0 for the relative probabilities of the final channels dW{s;a.az)- 
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